Chiral Symmetry Restoration and Scalar-Pseudoscalar partners in QCD 
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We present a consistent description of Scalar-Pseudoscalar partner degeneration at the QCD 
chiral transition in terms of the dominant low-energy physical states for the light quark sector. 
First, we obtain within model-independent one-loop Chiral Perturbation Theory (ChPT) that the 
QCD pseudoscalar susceptibility is proportional to the quark condensate. We also show that this 
chiral-restoring behaviour is compatible with recent lattice results for screening masses and gives rise 
to degeneration between the scalar and pseudoscalar susceptibilities (xs, Xp) around the transition 
point, consistently with an 0(4)-like current restoration pattern. This scenario is clearly confirmed 
by lattice data when we compare Xs{T) with the quark condensate, expected to scale as xp{T)- In 
addition, we show that saturating xs with the (t//o(500) broad resonance observed in pion scattering 
and including its finite temperature dependence, account for the peak structure of xsiT) in lattice 
data and the associated critical temperature. This is carried out within a unitarized ChPT scheme 
which generates the resonant state dynamically and is also consistent with partner degeneration. 

PACS numbers: ll.lO.Wx, ll.30.Rd, 12.39.Fe. 12.38.Gc. 



Chiral symmetry restoration SUviNf) x SUA{Nf) — > 
SUv{Nf) with Nf light quark flavours has been a mile- 
stone in our present understanding of the Quantum 
Chromodynamics (QCD) phase diagram and Hadronic 
Physics under extreme conditions of temperature T and 
baryon density, as those produced in Heavy-Ion and Nu- 
clear Matter experimental facilities such as RHIC, CERN 
(ALICE) and FAIR. Lattice simulations support that de- 
confinement and chiral restoration take place very close 
to one another in the phase diagram. In the physical 
case Nf = 2-1-1 (0 7^ 77i„ = rud = ruq <^ m^) and for 
vanishing baryon chemical potential, they point towards 
a smooth crossover transition at pseudocritical tempera- 
ture Tc ~ 145-165 MeV [Hll], the results being fairly con- 
sistent with the 0(4) universality class [3], which would 
hold for two light flavours in the chiral limit rUq = 0. 
The crossover nature of the transition means in particu- 
lar that there is no unique way to identify the transition 
point, the most efficient one in lattice being the scalar 
susceptibility peak position, rather than the vanishing 
point for the quark condensate {qq)T, the order parame- 



ter, which decreases asymptotically with T for niq ^ 0. 

The equivalence with the 0(4) -^ 0(3) breaking pat- 
tern, which in its simplest linear realization takes place 
through the cr-component of the 0(4) field (ct, tt") ac- 
quiring a thermal vacuum expectation value and mass 
vanishing at the transition in the chiral limit, led to early 
proposals of tt— cr meson degeneration ("chiral partners") 
at chiral restoration ^ . Degeneration in the vector-axial 
vector sector {p and a\ states) as a signature of chiral 
restoration has also been thoroughly studied [F'. Nowa- 
days, we know that the a state is well established as a 
TTTT scattering broad resonance for isospin and angular 
momentum / = J = 0, known as /o(500) ^, which is 
then difficult to accommodate as an asymptotically free 
state. Therefore, in order to study chiral partner degen- 
eration in the scalar-pseudoscalar sector, it is more ap- 
propriate to analyze the corresponding currents, which 
can be derived from a chiral effective lagrangian without 
introducing explicitly a particle-like a degree of freedom. 
The scalar and pseudoscalar susceptibilities in terms of 
the corresponding QCD SU{2) currents are given by: 
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where q = (u, d) is the quark field, the pseudoscalar 
current P"(x) = qj5T°-q{x), the Euclidean four- volume 
J^d^x = JI^ dr J d^x with /3 — 1/T and (•)t is a ther- 



mal average. For xp, Parity invariance of the QCD vac- 
uum {{P°')t = 0) and isospin symmetry have been used. 
Z[s,p\ is the QCD generating functional with scalar and 
pseudoscalar sources coupled to the massless lagrangian 



in the light sector as —s{x){qq){x) + ipa{x)P°'{x), so that 
Z[mq,0] is the QCD partition function. 

Thus, should the scalar and pseudoscalar currents be- 
come degenerate at chiral restoration, XriT) and Xs(T) 
would meet at that point. The same problem has been 
studied in nuclear matter at T = [7] where it has been 
found that to leading order in nuclear density the pseu- 
doscalar susceptibility decreases as the quark condensate. 
Therefore, since xs is expected to increase, as a measure 
of the fluctuations of the order parameter, at least up to 
the transition point, it seems plausible that they meet 
near the transition. In an ideal 0(4) pattern, the match- 
ing should take place at the maximum of xs- 

Since P" has the quantum numbers of the pion field 
7r°, its correlators are saturated by the pion state at 
low energies. Actually, from the low-energy theorems 
of current algebra (equivalent to the lowest order predic- 
tion of chiral lagrangians) one has P" ^ 2BoFtt'^ (from 
PCAC theorem) with Bq = M^/{2mq) and F,M the 
pion decay constant and mass respectively, so that xp '^ 
4_BgF^G'7r(p = 0) with Gj^^p) the pion propagator. This 
is equivalent to saturate the correlator inside the integral 
in Eq.([2| with the dominant pion state. At this order 
one can take the free Euclidean pion propagator (interac- 
tions are suppressed at low energies) G{p = 0) = 1/M^, 
so that using also the Gell-Mann-Oakes-Renner relation 
JVPF^ — ~mq{qq), we get xp ^ —{qq)/''nq, as a first in- 
dication of the relation between the pseudoscalar suscep- 
tibility and the quark condensate. The latter result can 
actually be obtained formally as a Ward Identity (WI) 
from the QCD lagrangian [5^, in connection with the defi- 
nition of the quark condensate for lattice Wilson fermions 
[S]. However, both sides of the identity suffer from QCD 
renormalization ambiguities, so that this WI is formally 
well-defined only for exact chiral symmetry [9l [10]. It 
is therefore interesting to study how this identity is re- 
alized within Chiral Perturbation Theory (ChPT) [TT] . 
which describes the low-energy chiral symmetry broken 
phase of QCD in a model-independent framework. The 
previous current-algebra results are actually valid only 
at very low energies, i.e., to leading order (LO) in the 



ChPT expansion in powers of a generic scale p, denoting 
pion momenta or temperature, relative, respectively, to 
A^ ^ 1 GeV and Tc- In particular, the LO prediction for 
Xp is temperature independent, so it is not obvious that 
it can be simply extrapolated as, say, (qq) — ?> {qq)T- Ac- 
tually, all the quantities involved are expected to change 
with temperature due to pion loop corrections, namely 
M,(r), F^{T) and (gg)(r). 

Similarly, from Eq.(fl]), one can relate xs with the 
propagator of a "cr- like state" such that it couples lin- 
early to the external scalar source s{x) in an explicit 
symmetry-breaking term Csb = 2BoFs{x)a{x). With- 
out further specification about its nature and its coupling 
to other physical states such as pious, one already gets 
Xs ^ 4:BqF'^G„{p = 0), suggesting a growing behaviour 
inversely proportional to M^ as the sigma state reduces 
its mass to become degenerate with the pion. 

The Next to Leading Order (NLO) corrections to 
Xp can be obtained systematically and in a model- 
independent way within ChPT, where one can also cal- 
culate the scalar susceptibility xs to a given order only 
in terms of pion degrees of freedom. The price to pay 
is that we expect to reproduce only the behaviour of 
XS,p{T) for low and moderate temperatures. However, 
since xp is dominated by pions, whose dynamics are well 
described through ChPT, we expect to obtain a reason- 
able qualitative description of its T behaviour, whereas 
standard ChPT misses the peak structure of xs near the 
transition. We note in turn that the LO for xs vanishes, 
unlike that of xp- The ChPT NLO result for xsiT) can 
be found in [T21[T3]. For xp{T) we consider the effective 
lagrangian £2+^4 + ■■ ■, where £2n — (^(p^"), including 
their dependence on the pseudoscalar source p"" as given 
in [11]. We follow similar steps as in [13l [14], now for 
the pseudoscalar correlator. The LO comes from £2 only 
and the NLO includes: a) G,r corrections to the LO, 
including one-loop £2 and tree level £4 terms, b) pion 
self- interactions 0{ttp'^ x ir^p^) in £2, c) crossed terms 
£2 = 0(p"7r) X £4 = C(/7r) and d) £4 = ©(pVtt^) 
terms. The final result is: 



X^'^^^iT) 



^Bl 



F^ 

M2 



1 

32^ 



,(4/.i 



h) 



2M2 



ffi(M,r) 



-0{F-'^) = - 



{qq^'^P^lT) 



+ OiF-^), 



(3) 



where subleading terms are labeled by their F^ depen- 
dence, the first term inside brackets is the LO current- 
algebra 0{F'^), F and M are the lagrangian pion mass 
and decay constant, related to the physical values Af^r — 
140_MeV, F^ ~ 93 MeV by 0{F-^) corrections HI] and 
I3, hi are renormalized scale-independent low-energy con- 



stants (LEC) of £4 [TT| such that the result in Eq.(l3|) is 
finite and independent of the low-energy renormalization 
scale. The thermal loop function is given by: 
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FIG. 1: Left: Comparison between the pseudoscalar screening 
mass ratio and A; ^ for the lattice data in ';16, (masses) and 
|17| (condensate) with the same lattice action and resolution 
and Tc — 196 MeV. Right: Scalar susceptibility versus Aj.s ~ 
Xp(T)/xp(0) from the data in [I] for which Tc ~ 155 MeV. In 
both graphs A;,^ = r{T)/r{Q) with r = {qq) — {2mq/ms){ss). 



Therefore, we have obtained the WI connecting (qq) 
and xp to NLO in model-independent ChPT, including 
finite-T effects. Furthermore, the niq dependence can- 
cels in Xp(T)/xp(0) = {qq){T)/{qq){0), where only me- 
son parameters show up. To this order, Xp(^)/Xf (0) = 
1 - 3gi{M,T)/{2F^) so that the LEG dependence also 
disappears. Note that hi comes from a contact term in 
£4 and is therefore another source of ambiguity in the 
NLO condensate [TT] . 

We draw two important conclusions from the previ- 
ous results. First, Xp{T) has a critical nature, scaling 
like the order parameter instead of the much softer be- 
haviour Xp{T) - 1/M2(T) with M.,,{T) the thermal pion 
mass, given at this order in [15], as a naive extension 
M — )■ Mt^{T) of the leading term would suggest. One can 
actually observe a clear signal of this critical behaviour in 
the lattice analysis of Euclidean correlators, which deter- 
mine their large-distance space-like screening mass M^'^ 
in different channels [TB]. From Eq.(l2]) we expect xp = 
Kp{0) ~ {Mp° '^)~^ with Afp° '^ the pole mass associated 
to Kp{p), the pseudoscalar correlator inside the integral 
in Eq.([2| in momentum space. Assuming a soft temper- 
ature behaviour for Mp°^''{T)/M'"'{T), which is plausi- 
ble below Tc, we can then explain the sudden increase 

this channel [16] since 
1/2 ^, 



of Mlf{T)/Mlf{0) observed 
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we expect that ratio to scale like 

[{qq)iO)/{qq){T)]^^'^. As we show in Fig|lJ(left panel) the 
correlation between these quantities is notorious, given 
the uncertainties involved. Data are taken from the 
same lattice group and under the same lattice conditions 
[IHITT]. Note that in lattice works {2mq/ms){ss) is sub- 
tracted from (qq) in order to avoid renormalization am- 
biguities. Estimating the T = condensates from NLO 
ChPT [5T] , this effect gives a 6% correction and, from the 
lattice values, it is about a 15% correction near Tc- Apart 
from the screening versus pole mass and the strange con- 
densate corrections, one should not forget about the typ- 
ical lattice uncertainties, like resolution, choice of action, 
staggered taste breaking and large pion masses [J [T^ . 
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FIG. 2: Scalar versus pseudoscalar susceptibilities in ChPT 
and in our unitarized description. We show for comparison 
the lattice data of Figfl] (right). 



The second conclusion is that the decrease of xp ^nd 
the increase of xs ^s they approach the critical point, 
leads to scalar-pseudoscalar susceptibility partner degen- 
eration, which in an ideal 0(4) pattern should take place 
near the xs peak. Once again, this behaviour is observed 
in lattice data. In Figll] (right panel) we plot the sub- 
tracted condensate, expected to scale as Xp(r)/xp(0), 
versus Xs{T)/xp{^), both from the lattice analysis in 
jlj. The T = values are taken from ChPT. The cur- 
rent degeneration is evident, not only at the critical point 
but also above it, where those two quantities remain very 
close to one another. 

We now turn to our analytical results based on effective 
theories, which are plotted in FigJ2] ChPT to NLO gives 
an increasing xs{T), matching xp{T) at Td — 0.9Tc, 
where {qq)^^^'^ {T^) = xf'^'^ {T^) = 0. Although the 
ChPT expansion is not expected to reproduce the criti- 
cal behaviour, this matching is a mo del- independent pre- 
diction. Near the chiral limit M^r <C T, where critical 
effects are meant to be enhanced, the degeneration point 
Td = Tc~ SAU/Att + 0{M^/Tc), approaching the chiral 
restoration temperature in that limit. 

Since the scalar susceptibility is dominated by the 
I = J = lightest state, which does not show up in 
the ChPT expansion, let us consider its unitarized ex- 
tension given by the Inverse Amplitude Method, which 
generates dynamically in SU{2) the /o(500) and p(770) 
resonances at p = as poles in the second Riemann sheet 
of the corresponding partial waves. More details and pole 
results at finite T can be found in [TSl [TU] . Saturating 
the scalar propagator with the /o(500) thermal state and 
assuming that its p = mass does not vary much with 
respect to the pole mass, we have 
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the model-independent ChPT predictions. This normal- 
ization compensates partly the difference between the 
p = Q and pole masses. Under this approximation, the 
full self-energy real part is Mg = M^ — T^/A, where 
{Mp — iTp/2)^ is the pole position in p"^ [TO] and the 
self-energy imaginary part vanishes at p = 0. Note that 
the width Tp for this state is not small compared to Mp. 

We plot Xs(^) i^^ Fig|2l The result agrees with stan- 
dard ChPT at low T and improves remarkably the be- 
haviour near the transition. It develops a maximum at 
Tc ~ 157 MeV, very much in agreement with the lattice 
data peak position. Mp is monotonically decreasing with 
T while Tp increases first due to thermal phase space but 
at some point mass reduction dominates and Tp decreases 
[TO] . Furthermore, approaching the chiral limit by tak- 
ing the Mtt = 10 MeV poles from [TH] gives a vanishing 
MsiT) at Tc ~ 118 MeV and hence a divergent Xs from 
Eq.([5]) at Tc- Thus, we get the expected Tc reduction and 
stronger xs growth [3] near the chiral limit. 

The quark condensate cannot be extracted directly 
from the unitarized susceptibility. However, we can ob- 
tain an approximate description by assuming that the 
relevant temperature and mass dependence, as far as 
the critical behaviour is concerned, comes from pion 
loop functions as 6{qq)^{T,M) = BoT'^g{T/M) and 
5xs = Blh{T/M), with 5f{T) = f{T)-f{Q). This T/M 
dependence holds actually to NLO ChPT, as in Eq.Q. 
Then, from Eq. (fTj), since 5xs — '^dS{qq)/dmq, we get 



g{x) ^ g{xo) + 



y3 



dy for x > xo, (6) 



with To — xqM <C M a suitable low-T scale below which 
we use directly NLO ChPT, which has a better analytic 
behaviour near T — 0. The h function is obtained from 
the T dependence of xs in Eq.([5]). The resulting (qq)^ is 
plotted in Fig(2]with Tq ~ 12 MeV [23. The critical be- 
haviour is again nicely improved compared to the ChPT 
curves and is in better agreement with lattice data in 
that region. In addition, we obtain once more a scalar- 
pseudoscalar matching near the xs peak and hence of 
chiral restoration. The corresponding (qq)^ near the chi- 
ral limit (M^ = 10 MeV) is much more abrupt, vanishing 
and matching xs ^t the same point, as expected. 

Recall that in the above unitarized analysis, we are not 
performing a fit to lattice points. We just use the same 
LEC which generate the T = physical /o(500), p states 
[19] and then provide our results for the susceptibility 
and condensate. Fitting the LEC within their uncertain- 
ties would improve the compatibility of our unitarized 
curves with lattice predictions below and near the tran- 
sition. We do not expect our effective theory analysis to 
reproduce the chiral restoration pattern above Tc- 

Summarizing, we have shown that a consistent scalar- 
pseudoscalar 0(4)-like current degeneration pattern at 
chiral symmetry restoration emerges from lattice simula- 
tions and analytical effective theory analysis. The pseu- 



doscalar susceptibility scales as the quark condensate and 
becomes degenerate with its "chiral partner" scalar sus- 
ceptibility close to the scalar transition peak. The lattice 
data supporting this picture are well accounted for by the 
dominant physical states: pions and the /o(5G0) scalar 
resonance generated in pion scattering at finite tempera- 
ture. Although we have restricted here, for simplicity, to 
the Nf = 2 case, the analysis can be extended to Nf — 3, 
where the role of other scalar states such as the ao(980) 
can also be studied [2D]. 
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For standard ChPT we use the same LEC values as in 
[121113] . Their infiuence is more important in xs, due to 
the vanishing of the LO, than in xf, (w)- 
We find very small numerical differences changing Tq be- 
tween 10-60 MeV. 



